Synergistic Polaron Formation in the Hubbard-Holstein Model at Small Doping 
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We study the effect of dynamical Holstein phonons on the physics of the Hubbard model at 
small doping using the dynamical cluster approximation on a 2 x 2 cluster. Non-local antiferro- 
magnetic correlations are found to significantly enhance the electron-phonon coupling, resulting in 
polaron formation for moderate coupling strengths. At finite doping, the electron-phonon coupling 
is found to strongly enhance the non-local spin correlations, indicating a synergistic interplay be- 
tween the electron-phonon coupling and antiferromagnetic correlations. Although it enhances the 
pairing interaction, the electron-phonon coupling is found to decrease the superconducting transition 
temperature, due to the reduction in the quasiparticle fraction. 
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Introduction. Phonons significantly modify the proper- 
ties of conventional metals, giving rise to phenomena such 
as superconductivity and charge-density wave formation. 
These effects are understood through conventional theo- 
ries. However, the effect of phonons on strongly corre- 
lated metals is less understood. This problem demands 
investigation given the experimental evidence of signifi- 
cant electron-phonon (EP) interactions in correlated ma- 
terials such as the cuprates and manganites 0, 0] . 

Theoretical interest is motivated by the suggestion that 
electronic correlations greatly enhance the effective EP 
coupling 0. In this scenario, the interaction of the holes 
with their antiferromagnetic (AF) background results in 
band narrowing which enhances the effective EP cou- 
pling, driving the system to the polaron regime even for 
moderate EP couplings |lJ- Lll- LjI ■ Quantum Monte Carlo 
(QMC) calculations 6], for the Hubbard-Holstein (HH) 
model predict that the EP coupling enhances d-wave 
pairing in the large coupling regime. Other calculations 
for the HH or closely related t-J models, including finite 
size 

SUM 

and dynamical mean field approximation 
(DMFA) calculations [lfj, [llj , suggest that phonons con- 
tribute strongly to spin and charge ordering, pairing, and 
phase separation. 

In this letter we emp loy the dynamical cluster approxi- 
mation (DCA) LU1 to study the two dimensional HH 
model at small doping. We investigate the role of the 
EP coupling on its one and two particle properties. The 
DOS, charge susceptibility and unscreened local moment 
indicate band narrowing and polaron formation at mod- 
erate EP couplings. Since this implies reduced mobility 
of the holes, antiferromagnetism is strongly enhanced at 
finite doping, but not at half filling. We find that the 
effective AF exchange as well as the pseudogap observed 
in both the single-particle spectra and the uniform spin 
susceptibility are nearly independent of the EP coupling. 
D-wave superconductivity is suppressed by the phonons 
due to the reduction of quasiparticle weight. 



The HH Hamiltonian reads 
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with a nearest-neighbor hopping t, an on-site Coulomb 
repulsion U, dispersionless optical phonons with fre- 
quency ujq and an on-site EP coupling g. The dimen- 
sionless EP coupling is defined as A = 2g 2 / \2MujqW) 
and represents the ratio of the single-electron lattice de- 
formation energy E p = g 2 / (2Mu>q) to half of the elec- 
tronic bandwidth W/2 = At. We present results for small 
phonon frequency ujq — 0.3t, U — 8t and values of A in 
the weak and intermediate coupling regime, A < 0.7. 

To study the Hamiltonian we employ the DCA, a 
cluster mean-field theory which maps the original lattice 
model onto a periodic cluster of size N c = L 2 embedded 
in a self-consistent host. Correlations up to a range L c are 
treated explicitly, while those at longer length scales are 
described at the mean-field level. With increasing cluster 
size, the DCA systematically interpolates between the 
single-site DMFA result and the exact result, while re- 
maining in the thermodynamic limit. Cluster mean field 
calculations on the simple Hubbard model are found to 
reproduce many of the features of the cuprates, including 
a Mott gap and strong AF correlations, d-wave supercon- 
ductivity and pseudogap behavior [l3| . 

We solve the cluster problem using a quantum Monte 
Carlo (QMC) algorithm 14] modified to perform the sum 
over both the discrete field used to decouple the Hub- 
bard repulsion, as well as the phonon field x. The space 
of configurations of the latter field is significantly larger 
than the former, so the HH code requires significantly 
more CPU time than required for the Hubbard model. 
Thus, present calculations are restricted to clusters of 
size N c — 4. The Maximum Entropy method is em- 
ployed to calculate the real frequency spectra. 
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FIG. 1: Single-particle DOS for several values of A. The pseu- 
dogap width is negligibly dependent on A. Inset: the elec- 
tronic kinetic energy, E^ in — ^2 ktT ^{k)nk >rj , versus T. N(0) 
and the kinetic energy gain are suppressed with increasing A. 



Results. The electronic DOS at 5% doping is plotted 
in Fig. ^ With increasing A the DOS at the Fermi en- 
ergy, iV(0), is strongly suppressed. The pseudogap in the 
DOS, as measured by the peak to peak distance, does not 
depend much on A. However, the peaks adjacent to the 
pseudogap are strongly suppressed by the EP coupling. 
The EP coupling also reduces the separation of the Mott 
peaks which is a measure of U e ff- The reduction depends 
on the phonon frequency (not shown) and for small ujq 
we find that it is much smaller than that expected in the 
antiadiabatic limit (U e ff — U — XW when ujq — > oo). 

The inset in Fig^shows that the electronic kinetic en- 
ergy gain at low temperatures decreases with increasing 
A. The quasiparticle fraction, which will be discussed 
below (see Fig. 0}, also falls with increasing A. These 
results, together with the suppression of -/V(O) and the 
peaks adjacent to the pseudogap, which in previous calcu- 
lations were associated with quasiparticles 0, are con- 
sistent with a significant decrease in the effective elec- 
tronic bandwidth and a reduction of the coherent part of 
the electronic Green's function due to EP coupling. 

The local dynamical charge susceptibility is plotted in 
Fig. ffia). It develops a narrow peak with increasing A, 
suggesting the development of a nearly localized band 
of polaron charge carriers. When dynamical phonons are 
considered there is no self-localization, making it difficult 
to define a critical A c for polaron formation. For A = 
the width of the peak is of order of t, while for A > 0.5 
the width of the peak becomes smaller than the phonon 
frequency ujq which we take to indicate the crossover to 
a polaron regime. Note that the peak found in the DCA 
(solid lines) is much narrower than that found by the 
DMFA (dashed lines) for the same parameters [lfij . This 
is consistent with the argument Q that the effective EP 
coupling is enhanced by non-local AF correlations, which 
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FIG. 2: (a) The local dynamic charge susceptibility for sev- 
eral values of A. The dashed (solid) lines are DMFA (DCA) 
results. A narrow band of low energy charge excitations de- 
velops as A increases which indicates polaron formation, (b) 
Static charge susceptibility (compressibility) versus T for the 
same values of A. The strong enhancement indicates the ten- 
dency to phase separation. 



are present in the DCA but not in the DMFA. Polaron 
formation is also associated with phase separation, since 
the electronic band is narrowed, i.e. i e /y is reduced, and 
the effective AF exchange J e f /, as we will show below, es- 
sentially remains unchanged for small uio/t. An increased 
Jeff/teff leads to phase separation since the system can 
gain more energy from spin exchange than it can from 
the kinetic energy of the carriers |l7j . In accord with 
this argument, the bulk charge susceptibility (or com- 
pressibility) , shown in Fig. EJb) , is strongly enhanced by 
the EP coupling (for larger values of A than those shown, 
it diverges at low T) indicating a charge ordering transi- 
tion (phase separation). 

The dynamic cluster spin susceptibility is plotted in 
Fig. Ola) for Q = (0, w). Due to remnant AF spin corre- 
lations, it has a magnon peak at an energy proportional 
to the effective exchange interaction J e ff- The fact that 
J e ff increases very slowly with A is expected from per- 
turbation theory in t [18| and exact diagonalization re- 
sults |l9l | when uiq is small. J e f / is independent of A 
when o>o — > since the lattice cannot respond during the 
virtual hopping between near neighbors which lowers the 
energy of the singlet relative to the triplet. 

The unscreened moment /i 2 = ((n-f — njj 2 } is plot- 
ted in Fig. |3fb). In the DCA calculation, increasing A 
greatly increases the temperature dependence of /x 2 , but 
at low temperatures, /i 2 is essentially independent of A. 
In the DMFA calculation, at the same temperatures the 
EP coupling is found to significantly reduce due to 
the tende ncy of the phonons to reduce the on-site cor- 
relations jljJl ■ In the DCA, the reduction of the effec- 
tive hopping amplitude due to the larger effective EP 
coupling (i.e. enhanced polaron formation) counteracts 
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FIG. 3: (a) The dynamic cluster spin susceptibility with Q = 
(tt, 0). The location of the peak is a measure of the effective 
spin exchange J e ff- (b) The unscreened moment (i 2 — ((n-f — 
n l) 2 ) versus T, calculated with the DCA (solid lines) and the 
DMFA (dashed lines). The EP coupling strongly increases 
the temperature-dependence of /j, 2 but the low-temperature 
value is unaffected, (c) The bulk spin susceptibility versus T. 



FIG. 4: (a) Tat versus filling for several values of A cal- 
culated with the DCA (solid lines) and the DMFA (dashed 
lines). When n — 1 , the DCA Tjv increases slightly, consis- 
tent with a small increase in J e / / shown in Fig.^Ja). At small 
doping, Tat is enhanced by the EP coupling in contrast to the 
DMFA result where Tjv is suppressed, (b) The Matsubara 
quasiparticle fraction Zq — 1/(1 — lmS(K,irT)/nT) versus T 
for several values of A. (c) Superconducting transition tem- 
perature versus A at 5% and 15% doping. 



the reduction of on-site correlations, thus leaving the low 
temperature unscreened moment unaffected. 

The bulk spin susceptibility is plotted in Fig. El (c). Un- 
like the charge susceptibility, the low temperature bulk 
spin susceptibility is essentially independent of A. When 
the bulk susceptibility is divided by the unscreened mo- 
ment /j, 2 , the four curves collapse upon each other (not 
shown), with small deviations consistent with the slight 
increase of J e // with A. This indicates that the temper- 
ature and A dependence of the bulk susceptibility, is due 
predominantly to its effect on the unscreened moment. 
The emergence of the pseudogap is associated with the 
suppression of the spin excitations and it is interesting 
to notice the weak A-dependence of both the low-T bulk 
spin susceptibility and the pseudogap in the DOS. This 
suggests that the pseudogap energy scale is J e ff- 

The AF transition temperature Tn at half filling (n = 
1) is weakly dependent on A as shown in Fig. 01 a). It in- 
creases only slightly, consistent with the increase in J e f f 
with A found in Fig. CJa). However, when the system is 
doped, T\r increases dramatically with A since polaron 
formation suppresses the mobility of the holes that dis- 
rupts the AF order. Hence, not only do the AF corre- 
lations enhance the EP coupling Q and polaron forma- 
tion 0, but the EP coupling also strongly enhances an- 
tiferromagnetism in the doped model. This establishes a 
synergistic cooperation between the EP coupling and AF 
correlations. Note that at half filling and small dopings, 
unlike DCA, the DMFA predicts a reduction of Tjy asso- 
ciated with the decrease in U e ff. However, like the DCA, 
at larger doping the AF order is enhanced in DMFA due 
to the decrease of charge carrier mobility, although the 



effect is not as pronounced. 

It is unclear a priori if the EP coupling will en- 
hance superconductivity. For U > 6t the EP vertex 
at small momentum transfers increases with increas- 
ing U, hence contributing to a d-wave pairing inter- 
action 0. However, the EP coupling also suppresses 
N(0) as well as the Matsubara quasiparticle fraction 
Z = 1/(1 - Im£(-ftT,7rT)/7rT). It is plotted in Fig.gfb) 
for different A for K = (0,tt). As T -> 0, Z Q -> Z where 
Z is the quasiparticle fraction defined by the derivative 
of the real-frequency self energy. At low temperatures 
the quasiparticle fraction is strongly reduced for large A. 
The superconducting transition temperature is plotted 
in Fig. 0{c) versus A at 5% and 15% doping. The lead- 
ing pairing instability remains even-frequency d-wave. T c 
falls with increasing A, indicating that the reduction of 
N(0) and the quasiparticle fraction outweighs the en- 
hanced interaction. Alternatively, several authors have 
proposed that pairing in the cuprates could be driven 
by kinetic energy gain ji^J, • The suppression of the 
kinetic energy seen in Fig. ^ is consistent with the sup- 
pression of T c in such scenarios. 

While our results are valid at small u>q, qualitative dif- 
ferences are expected when it increases. However, note 
that in real materials the phonon energy is even smaller 
than our value, e.g. in cuprates ujq w O.lOi, 0.15£ 0. 

Several experimental features in the cuprates can be 
addressed with our results. The enhancement of the EP 
coupling due to AF correlations can explain the large 
A w 1.5 seen in the photoemission spectra in the un- 
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derdoped cuprates Q], whereas density functional calcu- 
lations predict a moderate EP coupling A ~ 0.4 [22|. 
The enhancement of the charge susceptibility (Fig. |2J) 
shows that EP coupling is relevant for the char ge i nho- 
mogeneities observed in many high T c materials |23j . 

There is evidence in the literature that AF correlations 
can enhance EP coupling, and mixed results regarding 
polaron formation and pairing. Exact diagonalization 
calculations with both adiabatic and dynamic 0, 
phonons, finite-size QMC |6j and diagrammatic QMC p| 
all suggest that AF correlations strongly enhance the ef- 
fect of EP coupling. Those with dynamic Q phonons find 
evidence that the phonons can enhance the equal time 
AF correlations. However, these calculations have an ex- 
tremely truncated phonon Hilbert space that limits their 
ability to properly describe polaron formation. Perhaps 
as a consequence, these calculations suggest that Hol- 
stein phonons enhance d-wave superconductivity. Bipo- 
laron calculations [Isl 0] suggest enhanced formation 
of mobile intersite singlets for intermediate EP coupling. 
However, the weak A-dependence of the spin suscepti- 
bility, Fig. 01(c), excludes this scenario. Bipolaron cal- 
culations, like exact diagonalization at small doping, are 
restricted to one or two carriers, while the DCA addresses 
finite doping in the thermodynamic limit. In recent weak- 
coupling functional RG calculations, Holstein phonons 
were found to reduce the AF susceptibility, as well as 
pairing (but certain lower symmetry modes enhance pair- 
ing) [25j- However, these results neglect polaronic effects. 

Conclusion. We study the effect of Holstein phonons 
on the properties of the two-dimensional Hubbard model 
using the DCA and QMC. At small doping we find that 
iV(0) and the quasiparticle fraction are strongly sup- 
pressed with increasing EP coupling. This, together with 
a strongly increased charge susceptibility and the forma- 
tion of a narrow band of charge excitations in the dy- 
namic susceptibility, are taken as strong evidence for the 
formation of polarons. The evidence is far stronger in 
the DCA than the DMFA results, consistent with the 
argument that non-local spin correlations strongly en- 
hance the EP coupling |3j. The spin properties cap- 
tured by the DCA show that while the AF exchange 
J e f f increases slightly, the primary effect of the EP cou- 
pling is to increases the temperature dependence of the 
unscreened moment /j, 2 while essentially preserving its 
low-temperature value. Despite the reduction of U e ff 
the phonons do not suppress the low-T moment, since 
the polarons suppress the hybridization of the sites with 
their environment. The AF transition temperature T/v 
at half filling increases slightly, by an amount commen- 
surate with the increase in J e //. However, at finite dop- 
ing TV increases dramatically, due to the suppression of 
the itinerant holes mobility by polaron formation. Thus 
not only do the AF correlations increase the effective EP 
coupling but also the EP interaction increases the AF 
correlations. We find that both the magnitude of the 



pseudogap in the single particle DOS and the pseudogap 
signature in the bulk spin susceptibility are not affected 
significantly by EP coupling, suggesting that the pseu- 
dogap energy scale is J e ff- Finally, although the EP 
coupling is expected to increase the d-wave pairing inter- 
action |6j , we find that the superconducting T c is actually 
suppressed due to the reduction of the quasiparticle frac- 
tion Z and -/V(O). However it is possible that couplings 
to other phonon modes with lower symmetry, may con- 
tribute an interaction which will enhance the d-wave T c . 
This will be a subject of future studies. 

We thank G. Sawatzky, T. Devereaux, C. Honerkamp, 
R. Scalettar and H.-B. Schuttler for useful discussions. 
This research was supported by NSF DMR-03 12680, 
CMSN DOE DE-FG02-04ER46129, ONR N00014-05-1- 
0127 and NSF SCI-9619020 through resources provided 
by the San Diego Supercomputer Center. TM acknowl- 
edges the Center for Nanophase Materials Sciences, spon- 
sored by the Division of Scientific User Facilities, U.S. 
Department of Energy. BM acknowledges the UND Com- 
putational Research Center, supported by NSF EPS- 
0132289 and EPS-0447679, and the hospitality of the Pa- 
cific Institute of Theoretical Physics. 



[3. 
[4 
[5 
[6 
[7; 
[8 

[9 
[10 

[11 
[12 

[13 
[14 
[15 
[16 



li- 
tis 

[19 
[20 
[21 

[22 
[23 

[24 
[25 



A. Lanzara et ai, Nature (London) 412, 510 (2001); T. 

Cuk et al., Phys. Stat. Sol. (b) 242, 11 (2005). 

A. J. Millis et ai, Phys. Rev. Lett. 74, 5144 (1995); J. 

M. De Teresa et al, Nature (London) 386, 256 (1997); 

N. Mannella et ai, Nature (London) 438, 474 (2005). 

J. Zhong et al, Phys. Rev. Lett. 69, 1600 (1992). 

P. Prelovsek et ai, Phys. Rev. Lett. 96, 086402 (2006). 

A. Mishchenko et al, Phys. Rev. Lett. 93, 36402 (2004). 
Z.B. Huang et ai, Phys. Rev. B 68, 220507(R) (2003). 

B. J. Alder et ai, Phys. Rev. Lett. 79, 3022 (1997). 

J. Riera and A. Moreo, Phys. Rev. B 73, 014518 (2006). 
T. Sakai et ai, Phys. Rev. B 55, 8445 (1997). 

Sangiovanni et al, Phys. Rev. Lett. 94, 026401 (2005). 
Capone et al. Phys. Rev. Lett. 92, 106401 (2004). 
H. Hettler et ai, Phys. Rev. B 58, R7475 (1998); M. 
Hettler et al, Phys. Rev. B 61,12739 (2000). 
Maier et ai, Rev. Mod. Phys. 77, 1027 (2005). 



G. 
M 
M 
H. 
Th 



M. Jarrell et al, Phys. Rev. B 64, 195130 (2001). 
M. Jarrell et ai, Physics Reports 269 No.3, 133 (1996). 
At low T the DMFA DOS and dynamic charge suscepti- 
bility develop phonon features at ui ~ ujq not found with 
the DCA at the temperatures studied. 
V. J. Emery et al, Phys. Rev. Lett. 64, 475 (1990). 
A. Macridin et ai, Phys. Rev. B 69, 245111 (2004). 
W. Stephan et al, Physics Lett. A 227, 120-126 (1997). 
Th. Maier et al, Phys. Rev. Lett. 92, 027005 (2004). 
E. Arrigoni et ai, Phys. Rev. B 69, 214519 (2004); A. 
Bianconi et al, J. Phys.: Cond. Matter 12, 10655 (2000). 
O. Andersen et ai, J. Low Temp. Phys. 105, 285 (1996). 
J. M. Tranquada et al, Nature (London) 375, 561 (1995); 
J. E. Hoffman et ai, Science 295, 466 (2002). 
A. La Magna et ai, Phys. Rev. B 55, 14886 (1997). 
C. Honerkamp et al, preprint, cond-mat/0605161 



